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Abstract 

t> We introduce a spin model which exhibits the main properties of a Kerr medium to describe an 

00 

00 intensity dependent conpling between a two-level atom and the radiation field. We select a nnitary 

irreducible representation of the su(2) Lie algebra such that the number of excitations of the field is 

CN 

bounded from above. We analyze the behavior of both the atomic and the field quantum properties 

CO 



and its dependence on the maximal number of excitations. 
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I. INTRODUCTION 



The Jaynes-Cummings model (JCM) [I] is one of the simplest quantum systems describ- 
ing the interaction of matter with electromagnetic radiation employing a Hamiltonian of a 
two-level atom coupled to a single bosonic mode. Due to its importance in laser physics 
and quantum optics, this model has been the subject of many recent investigations, both 
theoretical and experimental [2] . It has been observed j3] that the temporal behavior of this 
system is very sensitive to the statistical properties of the radiation field in the initial state, 
revealing pure quantum features which have no classical counterpart, such as vacuum Rabi 
oscillations [4], collapse and revival phenomena [5] and squeezing of the radiation field [5J. 
These predictions have been verified experimentally [7]. 

The JCM has been generalized in different ways. In some cases, the interaction between 
the atom and the radiation field is no longer linear in the field variables, i.e. intensity de- 
pendent coupling has been considered |8]. Other investigations relate the JC Hamiltonian 
to u (1|1) and osp (2|2) superalgebras Most of the nonlinear generalizations of the JCM 
are made by using appropriate g-deformed oscillators [TUl [TT] . In this framework, the appli- 
cation of the corresponding quantum algebras has proved useful in obtaining exactly soluble 
models [121 [13]. Additional generalizations of the JCM can be found in |14j . 

Nonlinear phenomena are ubiquitous in quantum optics. One of the simplest phenomenon 
is the Kerr effect, which occurs when the refractive index of a medium varies with the number 
of excitations of the field. Quantum descriptions of optical fields propagating in a Kerr 
medium reveal a number of interesting features, such as photon antibunching, squeezing 
and the formation of Schrodinger cats [USEE], which have no classical analogues. The Kerr 
medium has been recently considered in the framework of g-deformed oscillators [18] and 
the Moyal phase space representation [19J. Some authors have considered a system where a 
two-level atom is surrounded by a Kerr medium using a special case of g-deformed oscillators 

EBH- 

In this paper, we use an SU(2) spin model, which exhibits the main properties of the 
Kerr medium [21] , to describe an intensity dependent coupling between a two-level atom and 
the radiation field in the framework of the JCM. The main feature of the model is that it 
is formulated on a finite dimensional Hilbert space, so that the number of excitations of the 
field is bounded from above. In the limit when the dimension of the Hilbert space becomes 
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infinite, the su (2) algebra of the spin operators contracts to the Heisenberg-Weyl algebra of 
boson operators and the model coincides with the usual JCM. 



II. A KERR HAMILTONIAN ARISING FROM AN 517(2) SPIN SYSTEM 

It is well known that the Heisenberg-Weyl algebra can be obtained by a contraction of 
the su(2) algebra [22J. Taking this limit as a motivation, let us start by considering the 
su (2) commutation relations 

[S+,S-] =2S 3 , [S 3 ,S ± ]=±S±, &fc = £i±z5 2 . (1) 
In a given unitary irreducible representation j, we introduce the operators 

' ~ ^ = 4^ n = S 3+3 , (2) 



V2j x'V 

which act on a (2j + l)-dimensional Hilbert space of a particle with spin j [2 lj . The su (2) 
algebra in this representation is thus written as 

[b,n}=b, [b\h] = -b\ [6,&t] = i__ ( 3 ) 

which allows the interpretation of b* (b) as creation (annihilation) operators for the quanta 
labeled by the number operator n. The oscillator-like Hamiltonian, in terms of the spin j 
operators, is 

H = ^(S + S_ + S-S + ) = Hlj (n + \ - |) . (4) 

The analogy of this expression with the usual optical Kerr Hamiltonian is evident, but the 
main difference is that the spin number excitation is bounded from above: n <2j. This is a 
consequence of the angular momentum condition —j<m< +j, once it is written in terms 
of the eigenvalues of the number operator n. Clearly the energy spectrum has the twofold 
degeneracy E (n) = E (2j — n). Note that although this Hamiltonian can be easily obtained 
in the framework of g-deformed oscillators, we use an algebraic model which possesses a 
physical interpretation of quantum deformation based on the fundamental SU (2) group. 
The standard harmonic oscillator limit is obtained when j — > oo, in which case b and V* turn 
out to be the usual bosonic operators acting on an infinite dimensional Hilbert space, i.e., 
the nonlinearity and the maximum number of excitations disappear. We observe that the 



3 



relation between the number operator and the corresponding creation-annihilation operators 
is 



n 



1 



(2j + 1; 



;tfb 



which leads to the usual relation n = b^b when j — > oo. 

The solution of Heisenberg's equation of motion for spin operators is 

S^t) = e^[cos(m)4 -sm(nt)S 2 ], 
S 2 (t) = e i ^[sin(f2t)S'i + cos(Qt)S 2 ], 

s 3 (t) = s 3 , 



(6) 



where Cl = — 4). Thus, the time evolution of S (t) is a rotation around the z axis, 
but the precession frequency depends on the excitation number operator. This result is in 
complete analogy with that obtained in the case of the Kerr medium [2UI [21] . 

In molecular physics, this model correspond to a diatomic molecule approximated by a 
Morse potential and the eigenstates correspond to the U(2) D SU{2) symmetry-adapted 
basis [251 I2S]. It was shown that this type of Hamiltonian describes the interaction of a 
collective atomic system with the off-resonant radiation field in a dispersive cavity [27\ . 



III. SU(2) COHERENT STATES 

Nonlinear coherent states (NCSs) have been discussed in different approaches and arriving 
to different states (2HI ED] ■ The most representative ones are constructed in the framework 
of g-deformed oscillators [31] , in which NCSs are eigenstates of the g-annihilation operator. 
A generalization of the SU(2) coherent states given by a multiphoton Holstein-Primakoff 
transformation is discussed in Ref. [32], and their physical consequences are studied in the 
framework of the standard JC model. Here we describe the SU(2) coherent states corre- 
sponding to the single photon case in the SU(2) spin model and use them to study the 
modified JC model to be defined in the next section. 

The group-theoretical coherent state (see ref. |28j) for the unitary irreducible representa- 
tion of the su{2) Lie algebra corresponding to Eq.(j2| is: 

|«) = e* 8 *-**- \j, -j) , (7) 
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where a = |e~^, < # < 7r, < (f) < 2tc and the ladder operators S± select the vacuum 
state |0) from the states \j, m) as usual: S- |0) = £_ \j, —j) = 0. The natural phase space is 
the sphere of radius j; the spin coherent states are thus represented by spots on the sphere. 
This dynamics leads to Schrddinger cat states on the sphere [21], i.e., a superposition of 
several spots located "far" from each other. In the harmonic limit, the sphere opens to the 
phase plane and the model coincides with the quantum harmonic oscillator or, through a 
renormalization, with the usual Kerr medium. 

Using the disentangling theorem for SU{2) operators, we can rewrite Eq.([7j) in the fol- 
lowing form 

io = (i + iei 2 )- j E ( 2j ) t ij+m) \j,m), (8) 

m =-j \j + m J 

where £ = tan (|) e~ 1 ^. Then, the photon number distribution of the field prepared in the 
state |£) is 

and the mean photon number is 

l£l 2 

(n) = 2j |g , 2 . (10) 

This result shows that (h) is bounded from above. It is convenient to rewrite Eq. ^ in the 
following form 

Vn{X)= j X" (1 - X) 2j " n , (11) 

which is a binomial distribution in n with parameter x = ^ [0, 1]. Evidently V n (x) 
converges towards the Poisson distribution if (n) remains fixed when 2j — > oo. One of the 
most interesting features of the SU(2) coherent states is that they exhibit squeezing which 
depends on 2j. 



IV. TWO-LEVEL ATOM SURROUNDED BY A KERR MEDIUM. 

We start with a brief review of the standard formulation of the JCM. In the dipole and 
rotating wave approximation, the JC Hamiltonian for a system of a single atom interacting 
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with a single mode is given by [TJ 

11 
7~Ljc = -fibo(a/a + aa)) + -hcoocr^ + h\(a + a + a'<7_), (12) 

where \ fouj(a^a+aa)) and |ftu;o0"3 are the well-known energy operators for the field and atom, 
respectively. Here u is the field mode frequency and Uq is the atomic transition frequency. 
The coupling between the atom and the radiation field is described by hX(a + a+a'a^), where 
A is a coupling constant. Besides a and a) are the usual bosonic creation and annihilation 
operators for photons in the mode, which obey the Heisenberg-Weyl algebra. a + and cr_ 
are the usual rising and lowering operators describing the fermionic two-level atom and 
is the atomic inversion operator, which follow the standard pseudo-spin algebra. 

Now, let us introduce the spin model discussed in Sec. II into the JC Hamiltonian to 
describe a two-level atom surrounded by a Kerr medium. This is achieved by replacing the 
usual bosonic creation and annihilation operators by the correspondent spin-j operators, i.e. 
we now consider 

Uj C = ^fwj(S + S- + SLS+) + ^hu cr 3 + -^(o+SL + S+&-). (13) 

Next, we will concentrate on studying the dynamics of the system. We first split the Hamil- 
tonian into two parts : Ho (the energy operator in the absence of interaction) and V (the 
coupling interaction), such that 

Ujc = n + v, (14) 

where 

"H = ^hu(S + S- + S-S+) + ^-hu a 3 , (15) 

4 J Z 
V = -^(a + S_ + S + a_). (16) 

In the interaction picture generated by T-Lq, the Hamiltonian of the system is 

, , -Hat a .Wo* 

Vj(t) = e'-^Ve-'-^. (17) 

It is clear that &± and S± become time-dependent by means of the Heisenberg equations of 
motion. Then the interaction picture Hamiltonian becomes 

V x (t) = _^(£ +e «W&_ + S+e'^a^), (18) 



where 



Q fl = oj -uj[l---^- 
J 2j 



(19) 



gives rise to the generalized detuning frequency, which depends on the excitation number 
operator. 

In order to solve the equation of motion in the interaction picture, we first observe that 
Vx describes processes where a photon in the mode is annihilated while the atom is excited, 
or vice versa. Then, at any time t, the state of the system \ip (£)) is expanded in terms of 
the states \n) <8> |+) and \n + 1) <g) ]— ), where \n) is an eigenstate of the excitation number 
operator and |+) and |— ) denote the atom in the excited and ground state, respectively. 
Thus 



2j 



w (0) = (*) > n ) 1+) + K (*) ' n + x > 



->]■ 



(20) 



n=0 



The equation of motion for the state of the system in the interaction picture, i.e. 
ih^~ t \if) (t)) = Vi(t) \ip(t)), gives the following simple coupled set of differential equations 
for the probability amplitudes 



whose general solution is 
a n (t) = \ a n (0) 



b n (t) = { b n (0) 



cos 



COS 



ia r , 



ib r , 



r t 



n 



+ ' ( ' - 2 - I' 

A A /(n + l)[l-^)e 



iClntl. 



(21) 
(22) 



r t 



■ ■ 

i^— sin . 

r„ V 2 



^(0)^/1 -gsin^ 



r t 



i— sin , 

r„ V 2 



ia n (0) W 1 



sin 



r 2 ™ v 2 



i^l fit 



(23) 



where a n (0) and b n (0) are determined from the initial conditions of the system and 



'fi2+4A 2 ( n + i) i_ 



2j 



(24) 



is the generalized Rabi frequency. 

This set of equations gives us the general solution of the problem. In order to calcu- 
late some physical quantities of interest, we still need to specify the initial photon number 
distribution of the field. 
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Exact solutions of generalized extensions of the JCM in the Schrodinger picture have been 
calculated in Ref. [Hj. Our modified algebra (Eq|3]) is also included in this generalization 
with the appropriate choice of the structure functions. Unlike this approach, here we use a 
unitary irreducible representation of the su(2) algebra to describe a Kerr medium, taking 
advantage of algebraic methods. In addition, we now explore the physical consequences of 
the model. 

V. QUANTUM PROPERTIES 

Since the JCM is a composite system, which includes the field and the atom sectors, the 
quantum properties can be described using the joint density operator 

P^(*) = llK*)XlM*)l. (25) 

where the subscript A and J 7 refers to the atom and the field contributions. 

The reduced density operators are then sufficient for calculating the averages of any 
dynamical variables that belong exclusively to one of the components. The reduced density 
operators are 

PA (*) = Trjr [p AT {t)\ , p T (t) = Tr A [p AT (t)} , (26) 

with which we are able to evaluate the expectation values of any atomic operator O a (e.g. 
population inversion operator 03) and of any field operator Oj (e.g. excitation number 
operator h) using the expressions 

(O a ) = Tr A [p A (t) A \ , {O t ) = Trjr \p F (t) O t \ . (27) 

The matrix density of the atom has dimension two and, unlike the usual JCM, the matrix 
density of the field has dimension (2j + 1). 

A. Atomic population inversion 

Many interesting atomic quantum effects have been observed in the context of the stan- 
dard JCM. Perhaps the most notable is the periodic transfer of population between the 
ground state and the excited state. These transitions constitute the so-called collapse and 
revival phenomena (CR). Physically, CR occur because an atom within a cavity undergoes 
reversible spontaneous emission, as it repeatedly emits and then reabsorbs radiation [33] . 
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The general expression for the atomic population inversion is 

(a 3 (t)) = Tr A [p A {t) a 3 ] = £ [k Wl" ~ K (t)\ 2 } , (28) 

n=0 



where we have used Eq. (27). We can see that the temporal evolution of (03 (t)) is essentially 
a result of a summation of probabilities at different Rabi frequencies. 

In order to compare the differences between the usual treatment and the spin model, let 
us consider the atom initially in the excited state and the initial photon number distribution 
described by the £77 (2) coherent states, i. e. \b n (0)| 2 = and \a n (0)| 2 = V n (x). Thus 




x n (i - xf- n 



cos (IV) 



(29) 



We observe that this solution converges towards the usual solution of JCM as 2j approaches 
infinity, i.e. when the su (2) algebra contracts to the Heisenberg-Weyl algebra. Simultane- 
ously, the photon number distribution of the field converges towards the Poisson distribution 
when 2jx remains fixed. 

It is well known that the system under consideration is sensitive to the statistical prop- 
erties of the electromagnetic field. In the g-deformed extensions of the JCM, it is usually 
considered that the field is initially prepared in a g-deformed coherent state (an eigenstate 
of the q -deformed annihilation operator). On the other hand, in Ref. [H] the field is pre- 
pared in the multipothon Holstein-Primakoff SU (2) coherent state in the framework of the 
standard JCM. In this work we consider both situations: an £77(2) spin model to describe 
an intensity dependent coupling, together with its associated coherent state as the initial 
photon number distribution. 

Numerical results of the atomic inversion in the exact resonant case (u = uq = A), with 
(n) = 20 photons on average and 2 7 — > 00, are shown in Fig. 1. It can be seen that the 
temporal evolution exhibits periodic collapses and revivals. This phenomenon also takes 
place in the spin model; however, in this case the structure is more complex due to its 
dependence on the maximum number of spin excitations 2j. Figures (2a), (2b) and (2c) 
show the temporal evolution of (03 (t)) for (n) = 20 photons on average and 2j = 1000, 100 
and 50, respectively. We observe in Fig. 2a that the sequence of CR is essentially the same 
as in the limiting case (Fig. 1), but the scaled time Xt^ needed for the largest revival of 
the initial value (a 3 (0)) depends on the maximum occupation 2j. An estimate of the scaled 
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At 

FIG. 1: Plot of the temporal evolution of the atomic inversion (a 3 (t)) in the standard 
JCM in the exact resonant case. The atom is initially in the excited state and the field is 
initially prepared in the standard coherent state with (n) = 20 photons on average. 



time (2j » 1) in this case gives 

\t R « - — . (30) 
y/ X 2 + (2 + 2 3X ) (1 - X) ~ Vx 2 + (1 + 2j X ) (1 - X) 

which agrees with the numerical results. It can be seen that when 2j decreases the structure 
starts to deteriorate (see Fig. (2b) and (2c)). Due to the symmetry of the photon number 
distribution, i.e. V n (x) = ^2j-n (1 — x)i periodic collapses and revivals also take place for 
any value of Xi except for the limiting case % — 1- I* 1 this case, which does not occur in 



the standard JCM, the leading term in the sum of Eq.(29) is that for which n = 2j and 



therefore (03 (t)) = 1, i.e. the atom remains in the excited state. 
B. Photon antibunching 

Both theoretically and experimentally, there is an interest in a variety of statistical prop- 
erties of the electromagnetic field, including the distributions of possible field energies, the 
photon number variance and the Mandel Q parameter (normalized second factorial moment) 

, ((An) 2 ) -(h) 

Q{t) = - \l } - (31) 

These statistical variables offer quantitative measures of how much the field differs from a 
classical field. In particular, the Mandel Q parameter vanishes for a Poissonian distribution. 
It provides information about the tendency of photons to arrive in bunches: when Q > the 
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FIG. 2: Plots of the temporal evolution of the atomic inversion (<5"3 (t)) in the spin model in 
the exact resonant case. The atom is initially in the excited state and the field is initially 
prepared in the SU(2) coherent state with (n) = 20 photons on average. Frames a), b) and 
c) correspond to 2j = 1000, 100 and 50 maximum number of excitations, respectively 



photons are bunched (super-Poisson), while for Q < the photons are antibunched (sub- 
Poisson) [33] . It is well know that a sub-Poissonian statistics is a signature of the quantum 
nature of the field. 



Averages appearing in Eq.(31) can be obtained easily from Eq.(27). From figure 3(a) we 
can see that the photon number distribution oscillates between sub-Poissonian and super- 
Poissonian statistics when 2j » 1. Figures 3(b) and (c) show that as 2j decreases, Q 
becomes negative and the photons are antibunched. From these results we infer the purely 
quantum regime of the spin model for the electromagnetic field. 
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FIG. 3: Plots of the temporal evolution of the Mandel Q (t) parameter in the spin model in 
the exact resonant case. The atom is initially in the excited state and the field is initially 
prepared in the SU(2) coherent state with (n) = 20 photons on average. Frames a), b) and 
c) correspond to 2j = 1000, 100 and 50 maximum number of excitations, respectively 



C. Squeezing 



To analyze the squeezing properties of the radiation field we introduce two hermitian 
quadrature operators: 



J r 



J II 



x = ^=, y = ^. (32) 

One of the consequences of the commutation relation for these operators is the uncertainty 



M 

3 



. When either of these variances is less than 



ft 2 



relation <(A£) 2 ) ((Ay) 2 ) > 
(uncertainty associated with the coherent field, including the vacuum), the state is said 
to be squeezed EE]- In the usual JCM, the cavity field surrounding a two-level atom 
initially excited, exhibits a time-varying pattern of squeezing both in the short time regime 



12 



and during the revivals [6j [37] . Thus the dynamics of the JCM leads to the squeezing of the 
radiation field, although the effect is rather weak. 

The variances of the quadrature operators can be expressed through the mean values of 
the spin operators 

1. 



<(A*) 2 ) 
<(Ay) 2 > 



1 

2j 
1 

2j 



Re <^J+J_ 
Re / J+J_ 



Re ( Ji > - Re ( J 



(33) 



Im { J+ 



In figure 4 we plot the temporal evolution of the uncertainties as a function of the scaled 
time. We observe in figure 4(a) that ((A£) 2 ) exhibits squeezing periodically in both cases, 
for 2j = 1000 (blue line) and 50 (red line) maximum number of excitations. In figure 4(b) 
it can be seen that ((A?/) 2 ) is larger than the minimum even at t — 0. On the other hand, 
figures 4(c) and (d) show that in the longer time behavior the uncertainties remain bounded 
in both cases (for 2j = 1000 and 50) . However, they always remain significantly larger than 
the minimum value and thus the field is no longer squeezed. 



VI. CONCLUSIONS 



In this work we have introduced a spin model which exhibits the main properties of a 
Kerr medium to describe an intensity dependent coupling between a two-level atom and 
the radiation field. The model is formulated in terms of spin operators acting on a finite 
dimensional Hilbert space, so that the number of excitations of the field is bounded from 
above. We have analyzed the behavior of both the atomic and the field quantum properties 
when the atom is initially in the excited state and the field is initially prepared in the SU (2) 
coherent state in the exact resonant case. 

It has been shown that the atomic population inversion exhibits periodic collapses and 
revivals when 2j 3> 1, while with decreasing 2j the structures start to deteriorate. When the 
mean photon number is maximal (x = 1), the atom remains in the excited state. As regards 
to the quantum properties of the field, we showed that the photon number distribution 
oscillates between sub-Poissonian and super-Poissonian statistics when 2j » 1, while as 
2j decreases Q becomes negative and the photons are antibunched. In the same fashion, we 
find squeezing only in the short time regime. We are currently exploring applications of a 
field-theoretical extension of these ideas [38] ■ 
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(c) (d) 

FIG. 4: Temporal evolution of the uncertainties as a function of the scaled time. Frames a) 

and b) show the short time regime of ((A£) 2 ) and ((Ay) 2 ), respectively. In the same 
order, frames c) and d) show the longer time behavior. The atom is initially in the excited 
state and the field is initially prepared in the SU(2) coherent state with (n) = 20 photons 
on average. The maximum number of excitations are 2j = 1000 (blue line) and 2j = 50 

(red line). 
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